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Confining configurations are introduced into the standard instan- 
f^ , ton vacuum model. This drastically improves theoretical properties of 

^^ I the vacuum: instanton size density d{p) stabilizes at /) ~ 0.2/m, all 

chiral effects are formulated in a gauge-invariant way and quarks are 
confined. An interesting interplay of chiral and confining dynamics is 
observed; for the realistic values of parameters the Georgi-Manohar 
picture emerges with chiral radius Rch ~ p ~ 0.2 fm much less than 
,^ . confining radius Re ~ hadron radius ~ Ifm. In the limit Rch *C Re 

the chiral mass Mchip) is unaffected by confinement and can be taken 
in the local limit M^hip = 0). 

Different types of effective chiral Lagrangians (ECL) are obtained, 
5^ I containing all or a part of gluon, quark and Nambu-Goldstone-meson 

fields. The ECL are manifestly gauge-invariant and in the limit of no 
gluon fields coincide with those found previously. 

The problem of the double role of the pion - as a Goldstone meson 
or as a qq system is briefly disscussed using confining ECL with quarks, 
mesons and gluons. 
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1 Introduction 

Effects of chiral symmetry breaking (CSB) and in particular dynamics of 
Nambu-Goldstone bosons liave been studied for a long time even before the 
discovery of QCD [1]. 

The effective chiral Lagrangian (ECL) has been introduced [2] before the 
notion of quarks appeared in physics. Nowadays ECL and chiral pertur- 
bation theory is a powerful method [3] and only the loop cut-off problem 
and numerous phenomenological coefficients of ECL remind that this is an 
effective theory, not yet deducible from the first principles of QCD. 

What is the relation of confinement to chiral physics? At first sight two 
phenomena are not directly connected. This is explicitly so in the instanton 
model of chiral vacuum [4,5], where instantons do not bear confinement [6], 
but create zero modes [7] which are necessary for the quark condensate [8] 
and a fortiori for CSB [9,10]. At the same time there is a coincidence of 
chiral and deconfining phase transition temperatures in lattice calculations 
[11], and some attempts exist [12] to explain the fact in the framework of the 
instanton model. 

At present there are two disconnected approaches in quark physics, which 
treat confinement or chiral dynamics. The last approach uses ECL and chiral 
perturbation theory, where the notion of confinement is absent [2,3]. More 
microscopic among models of this type is the instanton gas or liquid model 
[4,5]. This model yields microscopic chiral dynamics of quarks and produces 
ECL with coefficients which are calculable within the model. However quarks 
are deconfined and gauge invariance in the resulting ECL is violated. This 
of course is the price for keeping only instantons and disregarding confining 
degrees of freedom in the vacuum. The first approach is exemplified by 
the so-called relativistic potential model (RPM) [13], where chiral effects 
are disregarded fully and all mesons including Nambu-Goldstone bosons are 
treated as bound state of qq interacting via linear potential. 

Recently a more general approach has been used for the qq system (for a 
review see [14]), starting from QCD Lagrangian, where due to confinement a 
string was shown to appear for nonzero orbital momentum L, while at small 
L one has RPM. In this way in [15] RPM has been derived from QCD. 

But chiral as well as spin effects have been disregarded in [15]. 

In this paper we start to consider systematically chiral and confinement 
effects on the same footing. To this end we consider QCD vacuum as con- 



sisting of instantons A^^^ plus confining configurations B^. Instantons create 
zero modes, which after being mixed and shifted due interaction give rise to 
CSB, chiral quark mass and ECL, as it was shown in the instanton hquid 
model in [10,16]. 

Confining background, superimposed on instantons, plays double role. 
First, and very important, it strongly modifies instanton density d{p) at large 
instanton size p. This stems from the fact that charge renormalization at 
large distances in presence of -B^ strongly modifies [17]: logarithmic infrared 
divergence disappears, and the charge freezes at large distances [17] 

Air 

^siR) ^ , , R-2+m- (1-1) 

where m is the lowest excitation mass of the order of 1 GeV. As a consequence 
g^i^p) in the renormalized instanton action Smst = -^r^ does not grow at large 
p and the instanton density converges [18] 

Thus confinement stabilizes the instanton size p at some relatively small 
value p r^ po, 

Pq IGeV-^ ~ 0.2/m (1.3) 

m 

This value is slightly less than obtained in the instanton liquid model [5,10], 
where the size of instanton is fixed by instanton interaction at p ^ 0.3/m. 
The important point here is that now po depends on fundamental properties 
of the QCD vacuum rather than on model- dependent properties of instanton- 
instanton interaction. Second, confinement introduces important dynamical 
effects in calculation of all chiral effects, e.g. in calculation of the coefficients 
of ECL. In what follows we write a new ECL modified by confinement which 
is manifestly gauge invariant. Here one can see explicitly how two dynamics, 
chiral and confining, interplay and in particular one can check some provi- 
sions of Georgi and Manohar [19]. It appears that indeed the chiral scale 
Ac SB ~ p^^ and the scale of confinement, coinciding with average size of a 
hadron Aqcd, are different and one may disentangle in some cases two dy- 
namics. In particular the chiral mass in the limit Aqsb *^ Aqcd becomes a 
local operator dependent only on chiral dynamics. On larger scale the chiral 



massive quark interacts with antiquark via a string and this creates a new 
quark mass -constituent quark mass. 

Thus the latter contains three ingredients: (i) current (or Lagrangian) 
mass, (ii) chiral mass, (iii) string mass. In magnetic moments or spin- 
dependent forces exactly this constituent mass enters. 

Another interesting aspect of the interplay between chiral (instanton) 
and confining dynamics is the problem of double role of pion - as a Nambu- 
Goldstone particle or as a bound state of quark and antiquark. We obtain an 
expression for the qq Green's function which contains both chiral and confin- 
ing effects (i.e. the string due to nonperturbative gluons and pion exchanges) 
and is an extension of our earlier result [14,20] where only confinement was 
taken into account, and of earlier result of [16, 21] where only chiral effects 
have been considered. One can visualize there two possible type of poles, 
as suggested in [19] when one expands in pion field, and their position is 
separated (parametrically a Nambu-Goldstone pole at m = and qq pole at 
f^a ~ \f^)- We show explicitly that the second type of pole at rua ( a "quark 
model pole" ) exactly cancels in the expression for the total Green's function, 
and only Nambu-Goldstone pole survives, while new poles of unified dynam- 
ics appear heavier than nia. This is in agreement with the scenario suggested 
by Georgi and Manohar [19]. 

The paper is organized as follows. 

In section 2 we study the quark propagation in the instantonic vacuum 
in presence of confining background B^. Methods of multiple scattering 
theory in Ad are widely used and the framework of the Dyakonov and Petrov 
approach [10] is modified to include 5^. In this modification all expressions 
are manifestly gauge covariant in contrast to the original ones of [10]. 

Meanwhile we have justified in Appendix A the accuracy of the usually 
done approximation - keeping only zero mode in the quark Green's function 
in the field of a given instanton: we show that omitted terms are of the 
order mp <^ 1 as compared to the contribution of the zero mode. Another 
important estimate in Appendix B concerns the shift of the zero mode eigen- 
value due to the presence of confining field B^. It is shown that the shift 
5\ is < AOMeV and satisfies (5A • p <C 1, thus making all the method of [10] 
applicable also for the realistic inclusion of confining field. 

In section 3 we are applying another approach using a specific effective ac- 
tion suggested in [21-24] to calculate chiral effects for any number of fiavours. 
To this end a new simple derivation of the effective action is given in Ap- 



pendix C, which takes into account B^. In the same section 3 we consider 
interplay of chiral and confining effects in the example of the quark chiral 
mass and explicitly show that these effects can be disentangled when their 
ranges are much different. 

The case of two flavours is considered in Section 4. Bosonization method 
of [21] is slightly modifled by confinement. One obtains here the same "gap 
equation" for the chiral mass as for Nf = 1. 

In section 5 based on bosonization results of the previous section we 
deduce ECL for quarks, pions and B^ and also for pions and 5^ only (after 
integration over quarks). 

We also obtain here and in the Appendix D an expression for F^^, which is 
a modification of that of [10,21] due to confinement. In section 6 we consider 
the qq Green's function and discuss contribution to it from confinement and 
Nambu-Goldstone modes. Here the double face of pion is discussed and some 
remarks on the OZI rule and OZI-violating mechanisms are made. 

The last section 7 is devoted to summary and an outlook. 

2 Quark propagation in the instantonic vac- 
uum with the confining background. 

In this Section we shall extend the method of [10] to the case when the QCD 
vacuum contains in addition to the gas of instantons and antiinstantons also 
some confining configurations, using mostly the method and notations of ref. 
[10]. Our goal is to write the quark propagator and effective action in the 
one-flavour case in the form where the gauge invariance and conflnement are 
present explicitly. 

We start with the ansatz for the vacuum 

TV 

A,ix)=j:Ai\x)+B,ix) (2.1) 

where A^ = A^_|_ + A^-_ is the total number of instantons and antiinstantons in 
the 4-volume V, A^ is the fleld of (anti)instanton in the singular gauge and 
B^{x) is the background conflning fleld which ensures the observed string 
tension a through its fleld correlators < F^i,{x)Fx„{y) > etc. [25]. Note 
that the instanton gas does not provide nonzero string tension, therefore the 



known value of a = 0.2GeV'^ fixes the normalization of the background field 

We define as in [12,26] the total quark Green's function S, "free" Green's 
function Sq and the quark propagator S''-*-' in the i-th (anti)instanton field: 

S ={-iD{A)-imr\ 5» = (-iZ)(5 + A«) - im)-\ (2.2) 

Sq = {-iD{B)-im)-\ 

We also introduce a complete set of real eigenvalues A^*-* and eigenfunc- 
tions u!^^ on a given center i: 

- iD{B + AW)m« = X^^\'^\ l<i<N;0<n<oo (2.3) 

So that S''-*-' can be written as 

n Xn — im 

with the definition 

U = S^- ^W (2.5) 

one can write the exact equations for the total quark propagator S [26] 

S = S,-Y.Q,k, (2.6) 

i,k 
Qik = i^i^ik ~ ti^Q /_^Qjk (2.7) 

Eqs.(2.6)-2.7) are exact for the the given decomposition of the gauge field 
in (2.1). One can now make an approximation for tj using the fact that the 
sum over large values of n in (2.4) is close to the "free" Green's function 
5*0(3;, y), since A„ ~ Vv^, n ~ p^ at large n and p^ ^< -B^^^ >. Therefore 
one can approximate tj by a finite sum 

n=0 An' -im 

In what follows we shall often keep only the lowest term n = 0, Ag ~ 0, this 
approximation was exploited in a series of papers [10,23,24]: 

5«(x,y) = 5o(x,y) + %^^^^^ (2.9) 

Aq — i'm 



We study the accuracy of the approximation (2.9) in Appendix A. We show 
there that the omitted terms are of the order of 0{mp) as compared to lowest 
mode contribution. Another important topic concerns the shift of eigenval- 
ues A^*\ defined in (2.3-2.4), due to the background field B^. We argue in 
Appendix B, that this shift is insignicant for the realistic values of gluonic 
condensate and instanton radius p. In particular, (Aq )^ is of the order of 
(30 — AOMeVy and therefore does not spoil the approximation (2.9). 

Insertion of the separable form (2.8) into (2.7) yields the following solution 
for S: 

Six,y) = Soix,y)- ^ ««(x)(^m - A + V);^,fcwW+(l/) (2.10) 

i,k;n,m 

where we have defined the matrices \mn,ik, Vnm,ik 

Kn^,^fc = f u^:^+ {z){-tD{B) - tm)ul^\z)d*z (2.12) 



by definition Vnm, {i = k) = . 

The form (2.10) coincides with that used in [10], when one keeps in the 
sum only n = m = 0, Aq = 0. 

In the rest of this section we shall calculate the quark Green's function 
S averaging (2.10) over vacuum configurations B^ and (anti) instanton po- 
sitions and orientations. We shall follow here the direct approach of [10], 
which makes explicit the appearing of the chiral mass of the quark, while 
in the next sections we follow another approach [21-24], where the notion 
of the effective action is introduced from the beginning. The new element 
which we shall obtain in this section is the gauge covariant form of the quark 
propagator with confinement taken into account. 

For the latter and for the averaging over 5^ and A^^^ one must define a 
physical gauge-invariant quantity associated with S{x, y) (which itself is not 
gauge invariant). 

The simplest quantity is the Green's function of one light and one heavy 
quark in the limit when the heavy mass is infinitely large. 

Ghl{x, y) =< trrS{x, y)T<^{y, x) > w (2.13) 



Here tr is a trace over color and Lorentz indices, V = 75, 757^, 1, 7^ etc., and 

ry 
^{y,x) = Pexp ig B^{z)dz^ (2.14) 

J X 

The integral in (2.14) is taken along the straight line connecting x and y- 
this is the remnant of the heavy quark Green's function. 

The angular brackets in (2.13) denote averaging over fields 5^, A'^^ defined 
as follows. 

< 0(4), B,) >^„ ^ = / n -^^dn,dfi{B)0{A^ , B,) (2.15) 

where Ri is the position of (anti)instanton, Vti is its color orientation, and 
dn{B) is the standard integration measure for the field B^, the specific form 
of it is not needed for our purposes. 

The total gauge transformation for the field A^ 

A,{x) ^ U{x){A^{x) + -d^)U+{x) (2.16) 

can be conveniently split into a homogeneous one for A^*^ and inhomogeneous 
for B„ 






4) ^ f/(x)4)[/+(x) (2.17) 



B^{x) -^ U{x){B^{x) + -d^)U^{x) (2.18) 

From(2.1) and (2.17-2.18) one returns back to (2.16). The color orientation 
for A|j*^ can be made explicit using 

4'ix) = n,A^al (2.19) 

where A^^^ is the standard singular gauge form [7,27] 
- 2 p^ 

It is clear from (2.17) that under global gauge transformation U the con- 
stant matrix ^2^ is simply "rotated" 

VLi -^ UQi (2.21) 



We now turn to eigenfunction -u^*) . Under gauge transformation it transforms 
as 

u^\x)^Uix)u(^\x) (2.22) 

To make gauge dependence in -u^ explicit, we write it in the form 

M«(x) = $(x,i?«)fi,v^„(x - i?«) (2.23) 

where 

$(x,i?«) = Pexp{ig r B^dz^), (2.24) 

and (fn is the form of solution in the singular gauge, e.g. for the (anti)instanton 
zero-mode solution one has [7,28] 



ami 



^(^) = - , . / .,3/2 ^ (2-25) 

with f^^ = -7|£^°'"(±i) and a.m spin and color S'f/(2) indices, and + and - 
referring to the instanton and antiinstanton zero mode respectively. 
Since <l>(x, i?*-*-*) transforms as 

$(a;, i?«) ^ U{x)^{x, /?(')) f/+(i?«) (2.26) 

one can satisfy (2.22) imposing on r2j the transformation law consistent with 
(2.21) 

a -^ U{R^^)^i (2.27) 

It is clear now that Vnm,ik (2.12) is gauge invariant, while S{x,y) in (2.10) 
transforms in standard way: 

S{x,y)^U{x)S{x,y)U+{y) (2.28) 

and Ghl (2-13) is gauge invariant. We do now a drastic approximation as in 
[10] to keep in (2.10) only terms with n = m = 0, and the lowest eigenvalue 
Aq. In accordance with [10] we define 

/ 1 X ^Si^ f Dik{R^'\R^''\ni,nk),i,k of one tipe , . 

\rh + v'"' ^rh^\ P,fe(i?», i?«, fi„ fifc),^. A; of diff. types ^^'^^^ 



where i, k of one type means both i and k are instantons or both antiin- 
stantons and we have used notation irh = im — Xq. Following the same line 
of reasoning as in [10] we obtain the following equations for D, P where we 
suppress arguments for simplicity 



ik 



-^V,k^ - ^ I d^R^'Un~V,, X ^.^, D,k (2.30) 

tm im IV J tm 1 — tmo 

Several comments are in order. First, P^- and -D^ can be considered as 
probability amplitudes for a quark to travel from a center i to a center j 
with all possible centers being on its way. Both P^- and Dij are defined not 
to contain (infinitely many, in principle) returns to the centers i,j, and those 
should be accounted for separately (note that (2.29) accordingly is a precise 
definition only for paths without returns). 
This is done introducing the factor 

oo 

(1 - im5)~^ = J2 i^^^T (2-32) 



n=0 



where imS is the amplitude of returning to the center only once, 

5 = DuiR^'\R^'\ni,ni) (2.33) 

Second, averaging over all R^^\ i^j with j ^ i,k is assumed in Dik, Pik- This 
can be factorized out only in the limit N^ ^ oo and this limit is assumed 
everywhere below. 

It is convenient to factorize explicitly the dependence on Vt in Di^, Pn. 

D^k = vl„X^tUdlU^,),,,v,p^ (2.34) 

P^k = vX^intuCnv (2.35) 

v^k = vXjntuv:',,nv (2.36) 

From (2.34) one concludes that the 6 does not depend on i in (2.33) and 
is a gauge-invariant quantity {Nc -^ oo), if one exploits the fact that it is 



10 



diagonal in spin indices a, (3 {cf (2.25)) 

6 = ^Tr^{Du{R^'\R^'\n,,n,))= (2.37) 

1 

1 

The color trace appearing in the last equality in (2.37) signals the gauge 
invariance of 5. We are interested in the limit m — > 0, therefore we introduce 
the reduced quantities (J, /, v where dependence on ffi and y is made explicit, 

Ne ffi^ ' Ne ffi^ 

y^u^^VN^y^u. ,^ \ (2.38) 

JSei m[\ — imo) 

and obtain equations not containing m any more 

jik ^ _yik ^ I yij(iji{i)djk (2.39) 



The fundamental role in the chiral mass generation is played by the so- 
called "consistency equation" [10] which in our notations is 

2N 
SPcd- = :^ (2.40) 

Using (2.39) it can also be written as 

~e = \i^f^P^(^^>^^>^^) + (2-41) 

+ 0{^reHr{VVV)} 
The first term in (2.41) can be written explicitly: 

SpciV^kVki) = Spc I d^zd''z'd*R'^^^^+{z - R^'^)^{R^'^Z)iD^{z,R^^^ x(2.42) 
x^{z - i?«) ■ ^+{z' - RkMR^''\ z')tm{z\ R^'^Mz' - R{i)) 

11 



We note that rf*^, p'^ and V^'^ as well as rf*^, p'^, V"^ are gauge covariant 
and transform as 

^jik^yik^^ik^ _^ f/(i?»)(f^l^^^d^'=)f/+(i?(^)) (2.43) 

In contrast to that 6 and e are gauge invariant, as can seen in (2.37) and 
(2.42). 

The measurable physical quantity associated with f,dis the heavy- light 
Green's function Ghl and we use (2.10), (2.13), (2.29) to write it in the form 
similar to that of Eq.(40) from [10] (we remind the reader that D*^ and F'^^ 
are defined without returns to the centers i and k and therefore the returns, 
i.e. factors like jz^;^, should be inserted in (2.10)). 

< GHLix,y) >5_^w=< ^{y,x)TSoix,y)T > - 

N ,\ 5 . 



2VNc im 1 — imS 
X < r$(|/, x)r^{x - Ri)<l>{x, Ri)<l>{Ri, y)^{y - Ri) >b,r, 

' a2^ a2 (2.44) 



X < r$(y, x)T^{x - Ri)(!>{x, Ri)dik{Ri, Rk)(^{Rk, yW{y - Rk) >b,r,,r^, 

+ {inst. — * antiinst.) 



N ,2, 



X < r$(y, x)v>^{x - Ri)<i>{x, Ri)fik{Ri, Rk)HRk, y'My - Rk) >b,r,j,r, 

— {inst. -^ antiinst.) 

When the field 5^ is put equal to zero, one comes back to the situation 
studied in [10]. In this case it is useful to work in the momentum space and 
one obtains as in [10] 

S{p) ~< GMp) >= J+^ig^ (2.45) 

with 

eN 
M{p) = ^^pV'(p) (2.46) 

12 



and the equation for e ("consistency equation", which one can call also the 
"gap equation") obtains from (2.40) putting all $ = 1, 

N J (27r)4M2(p)+p2 ^ ■ '> 

This is the equation (36) of ref. [10] in the limit m ^ 0. 

The corresponding partition function can be written as [21] 

ZqcniNf = 1) = const J DijD^^exp[- J —^^^^{p){p - iM{p))^{p)] 

(2.48) 
Note that the resulting theory of massive quarks with the chiral mass M{p) 
is not gauge invariant, as is clearly seen in (2.48). This serious defect is 
cured in (2.44) which is fully gauge invariant, but the notion of the quark 
Green's function S{p) and the chiral mass can not be made gauge invariant 
in contrast to Ghl- 

3 Method of effective action -the case of one 
flavour 

In this section we adopt another procedure - the method of effective action 
suggested in [21-24]. One can prove [24] that this method yields the same 
Green's function as in the more direct but tedious method of [10], which we 
have used in the previous section. 

The proof given in [24] contains intermediate formulas which are diverging 
for m — > 0, i.e. the partition function Z ~ m'^'^-t' . Therefore we give a 
different derivation in Appendix C and check at each point the accuracy of 
physical approximations. In particular we show in Appendix C that there is 
no divergency for m — > and that leading term in the limit mp <^ 1 is the 
same as obtained in [24] when 5^ = 0. 

As shown in Appendix C the QCD partition function for quarks in the 
field (2.1) can be written as (in the limit m ^ 0) 

Nf N 

X < n n(2^^/ - / d'xi,pb4\x) / d'y4^^{y)ibi,f{y)) >^,,n, (3.1) 

13 



Here / refers to the flavour (total number of flavours is Nf) the meaning 
of 2imf is explained in Appendix C eq. C.8, and Uq (x) is defined in (2,23) 
and we did not specify the difference between instantonic and antiinstantonic 
zero modes, which should be kept in mind in (3.1). 

We now consider systematically the cases of Nf = 1 in this section and 
the case of Nf = 2 in the next section. 

The case of Nf = 1. 



In this case the averaging procedure over Ri,Vti factorizes in (3.1) and 
one can introduce a two-fermion vertex (to stress the analogy with [21] we 
keep notations of that paper) 



W/'^^^4^ / ™„W 



Y± = / dRidVti / ilj+{x)iDu^(;>{x)d^x / d^u^'' {y)iDi){y) 



dRi^+{x)tD^{x - i?,)i(1^75)$(a;, Ri, y)^+{y - R^)tD^{y)dydx (3.2) 

where ^{x, Ri, y) = $(x, Ri)^{Ri, y) 

Following [21] we introduce identically integrations over A+, A_, r+, r_ to 
obtain 

Dfi{B)D^ljDilj+ / -± ^— dV+ dV.expW (3.3) 

J — CO ZTT J —CO ZTT J —CO J — CO 

w = f d'^xij+ibi) + i\+{Y+ - r+) + iA_ (yi - r_) + iV+ in ^ + A^_ in ^ 

J V V 

Integrating over r+, r_ by the steepest descent method in the thermody- 

' V 

d\^ d\ 



namical limit N± ^ oo, V ^ oo, -^ = const., one obtains 



Zqcd = f Dfi{B)-±—D^D^+expW (3.4) 

J Zn Ztx 

W = N+\n — — -N+ + N^ In — ^ - N^ + f d^x^^iD^ + i\+Y+ + i\_Y_ 
iX+V iX^V J 

The integation over dX^dX^ can be also done by the steepest descent 
method when N± is large [21]. 

Before doing that we integrate over ip,ip~^; using notation 

Y± = iJ+^^^Yi ■ ^l; (3.5) 

14 



we obtain 



(iA-L (iA_ 



Zqcd = J Df,{B)^^exp{W' + W"), (3.6) 

W = N+ln—^-N+ + N_ln——-N_ (3.7) 

Trln(-Z)2 + F+F_A+A_), 

where Tr means trace over coordinates, color, and Lorentz indices. Integrat- 
ing now over dX^d\_ one finds A_|_,A_ in the extremum of {W + W") to 
satisfy equations 

^ = Tr{F+F_A_(-Z)2 + F+F_A+A_)-^} (3.9) 

^ = Tr{F+F_A+(-D2 + F+F_A+A_)-i} 

For A^+ = A^_ = y, one finds A+ = A_ = ^^ , where e satisfies the "gap 
equation" 

— = ) Tr(^, ^ — 3.10 

Since according to (3.2), (3.5), one has for Y±, 

y±= J dRiib^{x - Ri)^{x, Ri, y)^{y - Ri)ib (3.11) 

the r.h.s. of (3.10) and hence e is gauge invariant (only closed "contours" of 
$(x, -Rj, y) enter under the sign of Tr). 

In case when B^ = 0, $ = 1, ID ^> d, (3.11) reduces to 

y±{^,y) = J ^,e^'^^-'\P^±{p)? (3.12) 



and (3.10) becomes the same as (2.47) 

d^p M\p) AVN^ 



(27r)>2 + M2(p) A^ 
15 



(3.13) 



with M{p) given by (2.46). 

Now back to the case B^ ^ 0. We integrate in (3.4) over dX+dX^ first, 
replacing in W the X± by their extremal values (3.9) We obtain the effective 
action for quarks in the form (A'^^ = 1) 



Zqcd = const / DiJi{B)DipDip^exp I dxdytp {x)[iD-5{x—y)+iM{x,y)]^lj{y) 

(3.14) 
where the nonlocal mass operator is 

eN f ^ 

M{x,y) = ^^^ J dR^iD^ix - RiMx,R,,y)^+iy - R^iD (3.15) 

and e is to be defined from (3.10). 

Note, that the effective action in (3.14) is now gauge invariant in contrast 
to (2.48), since from (3.15) one deduces that under gauge transformations 
M{x, y) changes as 

M(x, y) -^ U{x)M{x, y)U^{y) (3.16) 

M(x,y) as given by (3.15) is what one may call the chiral mass operator. 
Inclusion of background field B^ makes it gauge covariant, but now M[x,y) 
is dependent on the confining forces, more explicitly, M{x, y) contains the 
field -B^ and therefore depends on the string which connects quark with 
antiquark (or with the string junction in baryon). 

To study this dependence more explicitly consider again the heavy-light 
Green's function, calculated with the help of (3.14): 

<GHL{x,y)>B= f DijDij+Dfi{B)ij{y)<l>{y,x)ij^{x)e^^''^+'^^^ = (3.17) 



= tr,my,x){tD + tM)~;) 

We have neglected additional qq pairs (quenched approximation or large N^ 
limit). Using the Feynman-Schwinger representation [14] (3.17 can be rewrit- 
ten as 

/•OO 

< Ghl{x, y) >B=< tr,{iD - iM) / dsDze-^t^,{y, x)$(|/, x) >b (3.18) 

JO 

where 

^.(a:,i/) = P(e-^(*'-^^)$.(x,i/)), K = J^ ^dX , ^F =\a,,F,, 

(3.19) 
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and the ordering operator P ensures the proper insertion of the mass operator 
M^ into the phase factor $z(a;, y) where the subscript z refers to the contour 
of $2 taken along the quark path. 

In (3.18) the field B^ enters M and the closed contour C formed by the 
straight line and the quark path between x and y, as shown in Fig.l, together 
with insertions of the mass operator. 

The dynamics in (3.18) is defined by the Wilson loop with mass insertions, 
shown in Fig. 1, which obeys the area law: 

< ^.(x, y)^{y, x) >=< Wm{C) >~ exp{-SM ■ (t) (3.20) 

where Sm is the minimal area of the contour C with mass insertions. Eq. 
(3.20) means that a string is formed between light and heavy quark trajec- 
tories and this string is influenced by the insertions of the mass operator. 

To understand better these insertions one may look at the quark propa- 
gator and expand it in powers of the mass operator: 

(ti) + iM)-i = {ib)-i - {ib)-ld%iM{u, v)d'v{ib)-i (3.21) 

The nonlocality of M{x,y) as can be seen in (3.15) is of the order of 
the p - average size of (anti)instantons (the integrand in (3.15) behaves as 
u-R.K I _"r.|4 at (|a; — i?i|, ly — _Rj| ^ p). At the same time the average size 
of a hadron in (3.18) - the width of the contour C in Fig. 1 - is of the order 
of the confinement radius Re ~ 0.5 fm for lowest states or larger than Re for 
excited states. 

Therefore we have to distinguish two cases. 
i)p > Re. In this case the nonlocality of M is strongly influenced and interre- 
lated with the dynamics of the string - one cannot separate effects of chiral 
symmetry braking (the chiral mass M) and confinement. 
ii)p <^ Re. In this case we can effectively replace in (3.15) ^{x,Ri,y) by 
$(x, y) and hence rewrite M as 

M{x,y) = $(a;,y)|M(p)e^^(^-^)^ (3.22) 

Furthermore we can neglect nonlocality of M[x, y) on the scale of {x—y) ~ Re 
and replace the operator M(x, y) by 



M{x, y) ^ 6{x - y) / d{x - y)M{x, y) = 6{x - y)M{0) (3.23) 
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Introduction of (3.23) into (3.14) yields now a gauge-invariant expression 
with constant chiral mass and eq.(3.18) assumes the form 

< Ghl{,x, y) >= tr{iD - iM{0)) f dsDze-^-'^^'^^h""^ (3.24) 

where the effects of CSB and confinement are separated. The first produced 
the chiral mass M(0) which enters instead of current mass, and the latter 
ensures the string dynamics between the now massive light quark and heavy 
antiquark. We conclude this Section with discussion of the quark condensate. 
From (3.14) we have 

< qq >M= itr < {-it) - iM)-^ >b (3.25) 

In the case B^ = (3.25) reduces to 

f d'^p M(p) 

where M{p) is the Fourier transform of (3.15) when $ = 1, and is the same as 
in (2.46). We note that the form (3.26) coincides with that given in [10,21]. 
When B^ ^ one should calculate the original expression (3.25) which 
can be rewritten with the help of (3.24) as 

< qq >^= - < / trMDze-^^Mx.x) >b (3.27) 

Jo 

where the contour C in ipc is a closed trajectory of the quark, to be integrated 
over in Dz. From (3.26) one can notice that the integral is effectively defined 
by the distances of the order of the radius of instanton p. When this radius 
is assumed to be much smaller than the confinement radius Re (the case ii) 
above ), then the effects of confinement are unimportant and give only a small 
correction to (3.26). Thus we have two different situations: in computing qq 
system (wave functions and masses) one can use the inequality p <^ R^ 
and keep only M(0) instead of M(p) in the first approximation in Green's 
function; in computing < gg > on the other hand one should keep dependence 
M{p) and can neglect in the first approximation effects of confinement. 



4 The case of two flavours 

In the general case of Nf flavours in the effective action in (3.1) there enters 
a 2Nf vertex [21,24] 

N 

Y± = {-ff fd^R,dn., f[ fd^xiPpDu^'^ I d^yu^'\y)iD^f{y) (4.1) 

Integrating over dQi for Nf = 2 and taking the limit Nc ^ oo one obtains 

Y± = f d^RdetJ±{R) (4.2) 

where 

(J^(m),„ = f d.xd.mhfix] , 



[J±iR))f9 = I dxdyij+ix)hlTl5)K{x,y,R)ijg{y) (4.3) 



and 

K{x, y, R) = ib^{x - R)^{x, R, y)^^{y - R)iD (4.4) 

The interaction (4.2) is reminiscent of the 'tHooft determinantal interaction 
[7], and can be also compared to the similar term deduced in [28]. 

The difference in our case, as also in the B^ = case of [21,28], is that our 
interaction is nonlocal with nonlocality of the order of p - size of instanton. 
In addition, in contrast to [21,24] the term (4.2) is gauge invariant and takes 
into account effects of confinement. 

The effective partition function similarly to (3.4) is given by (up to 
unessential factors and replacing iX± ^ g± , R ^ u 

Zqcd ^ dg+ dg^DiljDilj^expW2 

r ^ 

W2= {Y^ ijpD^pf{u) + g+detJ+{u) + g^detJ^{u)}du (4.5) 

—Nj^lngj^ — N^lng^ 

One can introduce as in [21] the 2x2 fiavour matrices r^^ , r^ = (t, i),a = 
1, ...4, and use the identity 

i^a)f9i^a)f'9' = ^'^^fg'^gf' (4-6) 
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One can exploit the Hubbard- Stratonovich tranformation and write Zqcd 
through additional functions La, Ra [21] to obtain 

Zqcd = J D^{B)dg+ J dg^ J D^D^+DLaDRaexpW2 (4.7) 

W2 = -N+lng+ - N_lng- + f d'^x[ippDipf+ 

+g+Ll{x) + g-Rl{x) + 2(7+L,(x) J^(x)+ 

+2g^Ra{x)r{x)] 

where 

Jl{x) = {r-)fg{Mx))fg (4.8) 

The effective action W2 contains both fermionic degrees of freedom ip, ip^ 
and bosonic La,Ra, and global parameters g^,g-. Integrating over DifjDip^ 
we get a fully bosonic effective action 

X = ln det ( '^ ^ 2^^^^ ] , (4.10) 

V 2g^RK iD J ' ^ ' 

where L = LaT^ , R = RaTa , and e.g.{LK)xy = J L{u)K{x,y,u)du. 

We shall now study following [21] the phenomenon of CSB using variables 
La,Ra- We shall look for the condensate of L^jR^, namely introducing 

o-(x) = L4{x) + Ri{x), (4.11) 

we show that 

< CT >^ 0, < L4 - i?4 >=< Li >=< i?i >= (4.12) 

Since L and R transform under SUl^Nj) x SUji{Nf) as 

L^Ul-LU^, R^UnRUt (4.13) 

the nonzero value of < o" > signals CSB. 

The insertion of (4.11) , (4.12) into (4.10) yields 

X = Tr ln{-D^ + g+g^a^K^) (4.14) 
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The integration over dg+dg-da via the steepest descent method yields 
equations for determination of (7_,_ , g_\, a*^"-* 

T = -'-±^^Va - 2ag,g_Tri . ^' ) = (4.15) 

dW iV± Va^ K' 

H a g^Tr{- 



dg± g± 4: ^ -D^ + g+g_a^K^ 

From (4.15) we find that go = —y^, and 

K'^ N 

aWoTri . ^ ) = +- (4.16) 

Comparing with the gap equation (3.10) we see that (Tq = ^^ (our nor- 
mahzation differs from that of [21], note also misprints in numerical coeffi- 
cients in Eqs. (27-29) of [21]). 

We can now again define the quark effective action if we insert in 14^2 in 
(4.7) the extremal values of La, Ra, g+, g- found above. 

We obtain as before the effective action (3.14) when we express o"o and 
go in terms of e. 

Thus the one-quark situation in case of Nf = 2 is the same as in the case 
Nf = 1: the effect of CSB is to create the chiral quark mass. 

This is true however only in the approximation (4.12) when all the effects 
of bosonization reduce to the creation of nonzero o"o, and no boson exchanges 
(quantum boson fields La, Ra) are allowed. In the next section we shall 
discuss these effects and derive effective chiral Lagrangian with confinement 
taken into account. 



5 Effective chiral Lagrangian and quarks 

We are now in position to calculate the "bosonic" effective Lagrangian W{L, R) 
in (4.9), inserting there extremal values of g+ = g- = go and < a >= ao, 

namely 

eN 2N, 
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To parametrize the quantum bosonic fields in L, R we introduce as in [21] 
the forms 

L = ^a,'-±^UV, R = ^ao'-^^VU^ (5.2) 

where U = expiiniTi) , V = expiiaiTi) ,i = 1,2, 3. The eight bosonic fields 
TTj, (Tj, cr, Tj correspond to eight field L^, i?a, a = 1, ...4. 
Insertion of (5.2) into (4.9) yields 

N r 
W{L, R) = ^J d'x{a\x) + rf{x))- (5.3) 

Tr ln{ib + i(l + a + r])UVM+ + i{l + a - r])VU+M.} 

where M± = M^^^ and M is defined in (3.15) and Tr is taken over coordi- 
nates, color and Lorentz indices. One should keep in mind that terms linear 
in fields should be suppressed since they vanish due to the steepest descent 
condition, yielding < r] >=< -Wi >=< ai >= 0. 

It is instructive to expand the effective action (5.3) in ttj, o"j, o", r^ and to 
find the corresponding quadratic terms yielding masses of mesons. 

For pions this procedure looks like 

-W{tt) = Tr ln{ib + iUM+ + iU+M^) = 
Tr Iniib + iM + A) = Tr Iniit) + iM) + Tr{{ib + ilViy^A)- (5.4) 
— Tr[{ib + iM)-^A{ib + iM^^A] = (-VT^ + W^'^'^ + W^^'^). 

where 

A = i{e'^^ - 1)M+ + iie-'"^^ - 1)M_, M = M+ + M_ (5.5) 

The first term, W^'^\ contains no pions and describes contribution of qq pairs 
interacting with background field B^, the quark being already massive due 
to appearence of chiral mass M. 

The second and third term, W^"^^ and W^^^ are depicted in Fig. 2 (a) and 
(b) respectively. 

Expanding in (5.4) A in powers of ttj and keeping only quadratic terms 
one obtains 

+ W^'^\rr) + W^^\7r) = (5.6) 
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where 

N(k, k') = -Tr\ . ^ . iMik + k') + . ^ . iMik) J- ^iMik')] 

2 'iD + iM iD + iM -iD + iM 



and 



< x\M{k)\y >=^—jK{x,y,u)e''''du ; M(0) = M . 

It is easy to see in (5.6) and we show it exphcitly in Appendix D that for 
7ia = const i.e. for k = k' = the sum in the square brackets on the r.h.s. of 
(5.6) vanishes, signalhng vanishing of the pion mass, as it should be for the 
spontaneous CSB with zero current quark masses. 

A similar analysis for other mesons can be done and results are similar 
to those of [21] . We note that masses of all mesons, other than ttj, do not 
vanish and are of the order of typical hadron mass. 

We turn now to the final topic of this section; effective chiral Lagrangian 
for pions. In this case we integrate out all meson degrees of freedom ex- 
ept that of the pion, since pion is the lightest particle dominating at small 
momenta. From (5.4) we obtain 

W{7r) = -Tr Iniib + iMU^) (5.7) 

where 

^22 ^ ' 

The expression (5.7) reduces to the chiral Lagrangian obtained in [21,16] 
when M -^ M(0) and 5^ -^ 0, so that it) -^ id. 

In the form given in (5.7) the gauge invariance is seen explicitly (for that 
the factor ^{x,R,y) should be kept in M as is given in (3.15)). 

The effective action (5.7) describes quark pairs (integrated out) propa- 
gating in the confining gluon field B^ and the chiral field f/5 - the remaining 
degrees of freedom - those of pions (ttj) and gluons (-B^). 

We can compare (5.7) with the standard chiral Lagrangian [2,3] 






J d'^xTr A^A^ + ^^ J d^xeapySsTrA^ApA^AsA, + ... (5.9) 



where 



A^ = U+td^U. 
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Comparing (5.6) and (5.7) with (5.9) we deduce the expression for the 
pion decay constant Fj,. To this end one should average in Zqcd (4.7) over 
fields -B^, which in the leading N^. order results in averaging over B^ the ef- 
fective action W{-n) (5.6-5.7). Since N{k, k') includes trace over coordinates, 
one has < N{k,k') >b= {27r)^5{k + k')N{k), and N{k) = k^No{k). Finally 
comparing (5.6) to (5.9) we have 

^F^ = Noik = 0) (5.10) 

When Bfj, is absent, we come back to Eq.(68) of [10]. For details see Appendix 
D. 



6 Correlation functions and the double na- 
ture of pion 

Our starting point now is the QCD partition function where we keep quark 
and meson degrees of fredom (in addition to B^). To obtain it , we integrate 
(4.7) over dg+dg-, which effectively reduces to insertion into W2 in (4.7) 
g+ = g- = go- We obtain 

W{L,R,iIj)= fd^xd^x''^j{x)[iD + i(l + a + 7])UV^—^M+ (6.1) 

+ {l + a- .^)\/f/+i±^M];,,..,^,(x') + ^J d'x{a\x) + rf{x)) 

where t/, V are given in (5.2). 

We now can calculate the qq correlation function 

Tf{x,y) =< ■^+{x)rilj{x)'^+{y)rilj{y) > (6.2) 

where F = 1,75,7^)7^75, o-^^.- The result is 

n^(a:, y) = ^J D<l>,Dfi{B)e-'^^'^''''>{Tr{S{x, y; <D, B)TS{y, x, $, B)T)- 

(6.3) 
Tr S{x, x; $, B)Tr S{y, y; $, B)} 
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Here W{L,R) is the effective action (5.3); $j(a;) i = 1,2, ...8 denotes the 
set of meson variables; a{x),r]{x),aa{x),7ia{x) . S{x,y;^,B) is the quark 
propagator in the external field of mesons {$} and confining gluons {B^}. 

S{x, y; $, B) = [iD + i{l + a + r])UV^—^M+ (6.4) 

+^(1 + a- r^)\/[/+i±^M]-^ 

Two contributions in eq. (6.3) are shown in Fig.3a,b, we shall call them one 
(quark) loop and two-(quark) loop respectively. This however is somewhat 
misleading since the action W{L, R) corresponds to the superposition of any 
number of quark loops , since 

TV r 
W{L, R) = -Tr InS-^ + -f? / ^^^^(^^(a;) + 7f{x)) (6.5) 

Therefore our terminology refers only to the interior of the curly brackets 
in (6.3) with understanding that these graphs should be integrated with the 
weight shown in (6.3). Now we discuss the result of integration of diagrams 
(a) and (b) of Fig. 3 over bosonic fields in (6.3). To this end we expand 
W{L, R) = Vr($i) in powers of <l>j up to the quadratic terms [21]. 

W{^i) = - j dxdy ^^{x)Ci{x, y)<^,{y) (6.6) 

The term £j gets contribution from two diagrams depicted in Fig. 2(a,b). 
For the case of pions C^^ can be read off from (5.6). For small values of 
momenta the Fourner transform Ci{k = 0) plays the role of meson mass. As 
we discussed in section 5 contributions to CT^^k = 0) of Fig. 2(a) and (b) 
cancel each other establishing in this way the Goldstone theorem. One can 
easily see that Ci{x,y) do not depend on Nc (since also M does not depend 
on Nc). On the other hand. 



N g 

V 



2 



57r 



< F' >,„,,- 0{N,) (6.7) 



Since W^($) ~ ^ ~ N^, the exchange of a meson yields correction 0(1/A'"c 



N 
V 

- in other words one can tell that the coupling constant of a quark to the 
meson 
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Hence the contribution of the diagram Fig. 3(b) is 0(A^^ ■ ■^) = 0{Nc) just as 
that of the diagram Fig. 3(a). This is in contrast to the situation with gluon 
exchanges for the diagram of Fig.4, which contributes 0(A'"^ ■ ■^) = 0{N^) - 

^^ c 

this is the so-called OZI supressed diagram, just as the diagram of Fig. 3(b). 
However in the vector and tensor channels the latter does not exist (tensor 
mesons appear in the 1/Nc corrections, see [21,24]) and therefore in these 
channels OZI suppression can be explained by the smallness of the diagram 
Fig.4. At the same time, in the scalar and pseudoscalar channels there is no 
suppression since the diagram of Fig. 3(b) provides for the amplitude roughly 
the same amount as the OZI allowed diagram of Fig. 3(a). This conclusion 
of no OZI suppression in scalar and pseudoscalar channels has been obtained 
before in [29] using another arguments. 

We now turn to the important question about the double role of the pion 
- as the Goldstone particle and as the qq bound state. Our present formalism 
is a convenient setting for the study of this question and we shall actually 
have in mind all scalar and pseudoscalar particles (to be compared later with 
vector and tensor particles). 

Let us start with the effective action (6.1). It contains both quark and 
meson degrees of freedom, but here they enter not on equal footing. Quarks 
ipfyipf are supplied with kinetic term and are respectable quantum field- 
theoretical quantities. In contrast to that meson fields <l>j enter as auxiliary 
fields, they do not possess kinetic terms and should be considered as external 
fields to be averaged out with the given weight. When one integrates out 
quarks, one obtains an effective chiral Lagrangian as in (5.7-5.8), where the 
chiral pion field now inherits the kinetic term and the full QFT status. But 
in doing that (integrating out quarks) we average over all quark structure of 
the pion and "see only its chiral face". Quantitatively both faces of the pion 
are given by (6.3), or diagrams Fig.3(a,b). While the diagram of Fig.3(a) is 
the usual one considered in quark models, the diagram Fig. 3(b) is pertinent 
to the chiral degrees of freedom - it contains the chiral propagator {Ci)~y. 

Now in the limit of large N^ the diagram Fig. 3(a) will contain only gluon 
(B^) exchanges and provides only poles at p^ = m'^{n),n = 1, 2.... The chiral 
propagator {Ci)~y provides a pole at p"^ = ml; for pion mf = 0. How these 
two types of poles coexist? 

For pion when p^ is small, p"^ <^ Rq'^j P~^ the situation is relatively simple. 
In this region only the diagram Fig. 3(b) has the pole and it is dominant. 
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When p^ increases and is close to the pole m^, the same pole appears in 
Ci{p). This is the signal that the expansion in powers of $j done in (6.6) 
and considering the diagrams of Fig. 3(b) is meaningless - the notion of the 
meson mass is not useful when it is strongly p - dependent and even has a 
pole at p^ = ml. Formally, the diagrams Fig. 3(b) acquire a pole at ml, so 
that the sum of contributions of Fig.3a and Fig. 3b according to (6.3) can be 
written as 

Mk) = nf\k) - 4'\k)^^^i'\k) (6.9) 

where tt^ (k) is the contribution of the quark loop of Fig. 3a. 
Now N{k) can be written as (in the limit p ^ 0) 

N{k)=C + M^{0)4^\k) (6.10) 

The "quark model poles" appear in n^ '{k) , 

^f\k) = -^^ (6.11) 

k'^ + m^ 

One can see in (6.9) (at least in the limit p — *> 0, i.e. for kp <^ 1) that "quark 
model poles" (6.11) exactly cancel, and only Nambu-Goldstone pole at /c = 
due to A^(A; ^ 0) = |F^ ■ k"^ survives. Thus the Goldstone theorem manifests 
itself in our case, as it should be since chiral symmetry breaks spontaneously. 

7 Conclusions and outlook 

We have shown that the QCD vacuum with confinement and topological 
charges can be adequately described by the model, in which confining con- 
figurations are added to instantons. 

On one hand, instantons are stabilized by confinement, and become more 
dilute, so that they can be properly treated in the instanton gas approxima- 
tion. 

On an other hand, quarks and gluons are confined in the model, and 
unphysical features of the instanton gas or liquid model, where quarks prop- 
agate freely, are now absent. 

Therefore in this paper we have obtained a realistic approach, in which 
quarks, gluons and Goldstone bosons can be treated simultaneously on the 
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fundamental level. The basic effective action is given in (5.3) and can be used 
for diagrammatic expansion in powers of Goldstone exchanges or integrating 
over all bosonic fields, in particular when there is an extremum corresponding 
to a selfconsistent bosonic field of solitonic type. 

The latter might be important for baryons. 

In the limit when confining configurations vanish we come back to the 
effective chiral Lagrangian obtained earlier [21], [30]. 

The framework suggested in the paper can be used for the calculations 
of all effects where chiral physics and confinement are both important. A 
systematic quantitative study is planned in subsequent publications. 
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Figure captions 

Fig. 1. Graphical representation of tlie lieavy-liglit Green's function, 
eq.(3.17), Tlie straiglit line from x to y corresponds to the heavy-quark 
path, the light-quark path is typically away at a distance of Re ~ Ifm from 
the heavy quark. 

The chiral mass insertion (second term in Eq. (3.21)) is shown near the 
instanton position Ri with nonlocality of the order of p ~ 0.2 fm. 

Fig. 2. Graphical representation of two terms in the effective action (5.4). 
Part (a) corresponds to W^"^^ and part (b) to W^^\ solid lines denote quark 
propagator {iD + iM)~^, broken lines - emitted pion field ttj from expansion 
of A, Eq. (5.5). 

Fig. 3. Graphical representation of the qq correlator Il^{x,y), Eq. (6.3). 
Fig. 3a corresponds to the "one-loop term" (first term inside the curly 
brackets of (6.3)), while Fig. 3b corresponds to the "two-loop term" (second 
term inside the curly brackets). Broken line denotes the boson propagator 
which appears when one expands S{x, x, $, B) in powers of $. 

Fig. 4. The OZI violating two-gluon exchange diagram obtained from 
expansion in 5^ of the two-loop term of Eq. (6.3). 
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Appendix A 

Quark Green's function in the instanton field 

In this appendix we remind the reader the expansion of the quark Green's 
function g in the pure instanton field A^''\x) in the unitary gauge, which was 
found in [31]. 

G{x, y) = {ib - im) A^^^ + i^D ■ l^Jl + — m(x, y) (A.l) 

^ ^ Lilt 

where u{x,y) is the contribution of zero modes 

u{x,y) =uo{x)u^{y) (A.2) 

while A{x,y) is the Green's function of scalar particles in the instantonic 
field in the unitary gauge 

{-Dl + m')A{x,y)=6'^'\x-y) (A.3) 

One can obtain expansion of A in powers of m [32,29] in the singular gauge 
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A{x,y)=<p-'^\x)A{x,y)r'^\y), 0=1 + ^ (A.4) 



Consequently one has for G{x, y), x ^^ y ^^ p 

G(x,y) = — u(x,y)(l + 0(mp) (A.6) 

This can be compared with the spectral decomposition (2.4), with the result 
that the contribution of the nonzero modes is finite for m ^ 0. 

Now we can see the physical parameter of expansion in the ansatz (2.9) 
used here and earlier papers [10,21-24]. 

One can state that the omitted terms in (2.9) are of the order mp (since 
mx and my which can also appear effectively enter our expressions for quark 
propagator or effective action at distances x,y < p). 

Thus mp <^ 1 is also accuracy of our approximation in (3.1). One can 
see that for a typical instanton size p ~ 0.2 fm, u, d and s quarks satisfy 
condition mp <^ 1, while for c quark this is already violated and all effective 
actions considered below are not applicable for the c quark. 
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Appendix B 

Shift of eigenvalues due to background B^ 

We study eigenvalues in the field of one instanton A^ plus background 
Bfj,, which are to be found from the equation [12]: 

[-{d^ - tg{A^ + B^))^ - ga{E{A + B) + B{A + B))]^^ = A^v9„ (B.l) 

where A^ = -rja^u ^fl 2 f? while E,B are colorelectric and colormagnetic 
fields respectively. 

Let us concentrate on the shift of the zero eigenvalue which we shall 
evaluate by perturbation theory 

6X^ = (^o^V^^o) (B.2) 

where 6V is to be read off from eq. (B.l) and is due to B^ 7^ 0. One can 
extract the phase factor out (/?„, as is done in (2.23) and the rest is equivalent 
to the wave function in the Fock-Schwinger gauge for 5^; 

B^{x) = / a{u)dUyF^^{u) , a{u) = - (B.3) 

JO X 

In SV there are three terms: (i) linear in B^ (ii) quadratic in 5^ , g^B^^ 
(iii) proportional to E + B. Due to symmetry reasons only quadratic term 
contributes to SX^. 
We have 



/o (p^ + r2)3 

^00 y'^dv fpy 
= 2g^ / -— — r-- / a{u)a{u)duydu^ < F^^{u)F^,^{u) > 
Jo (1 + y'^y^ Jo 

Here we have introduced for an estimate the average value of < F{u)F{u') > 
which at -u = -u' should be less or equal to the gluonic condensate of [33]: 

^ < f;,f;, >= o.ouGeV (b.s) 

Taking into account that < F{u)F{u') > falls off at distances \u — u'\ ~ 
Tg ^ 0.2 fm [34], one can compute the integrals in (B.4) to obtain 



j2 
SX' ^ '—— < F^F" >< O.OOlSGet;^ (B.6) 
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In getting upper bound on the r.h.s. of (B.6) we use the fact that part of 
gluonic condensate should be due to instantons and in our case only that 
part of < F°'F"' > enters in (B.6) which is due to confining background B^. 
Hence we have average shift of the zero eigenvalue < \SX\ >< 40MeV. This 
is small as compared to p^^ ^ IGev , and therefore |(5A|p ^ 1. 
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Appendix C 

Derivation of the effective action (3.1) 

We start with the general form of the partion function for the quarks of 
flavour / = 1, ...Nf in the field (2.1) 

Z = const I Dfi{B)DijDip^exp - f ijjS'^^fdx (C.l) 

where the action of gluonic field is included in Dfi{B). 

We now transform S~^ in (C.l) to make explicit contribution of zero 
modes 

N 
i=l 

Using (2.5) and (2.8) after some algebra one obtains 

N K 

S-\x,y) = S,\x)5ix-y)+Y: E S,\x)ul:\x)itm-\-V)T,^^^u^^^^y)S,\y) 

i=l n,q=0 

(C.3) 
where matrix elements of V are 

Vn,,^^ = j uf^ {z){-lb{B) - im)uf {z)dz (C.4) 

It is easy to check that inverion of S~^ given by (C.3) yields (2.10). 

For n = g = Uq has a definite chirality and therefore Voo,m = —im. 
Hence, keeping only zero models in (C.3), as it is done e.g. in [21-24], one 
obtains a term 0( ^^^_^ ) in S~^ in (C.3) which was not present from the 
beginning in (C.2) 

The reason for this apparent paradox can be identified as an improper 
omission of terms with n, g 7^ 0. Indeed, matrix elements Vu^a and Viq^u are 
of the order of O(-), since eigensolutions of (2.3) with n 7^ have no definite 
chirality and have nonzero matrix elements even in the limit m = Q. 

Therefore {im—\—V)~^ ~ {det{im—\ — V))~^ is finite for m -^ 0, Aq ^ 0, 
and this fact resolves the apparent paradox occurring in derivation in [21, 
24]. 
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Insertion of (C.3) into (C.l) yields 

Z = const / D ii[B)Dil)D%l)^ exp{— / ^ljJSQ^^ljfdx)exp / Cdxdy (C.5) 
where 

N 

C = -J2i,jix)J2T.So\^H\^)i^rn - X - VU,,u^^^y)S,\y)tf{y) 

f 4=1 n,q 

(C.6) 
One can notice that / Cdx dy contains fermionic operators of the type: 

^j{x)S,\x)u';i\x)dx ^ ^^if) (C.7) 

Due to anticommutativity of ■^/'t the product ^}^*-'(/)\l/Jj*'*(/) vanishes. 
Therefore in expansion of expC only finite number of terms survives. Namely, 
if we for simplicity keep only zero modes, n = g = in (C.6) then one can 
write. 

N,Nf 

expjcdxdy = n (1 - ^?(/)(^^ - Ao - V^)-^^?+(/)) (C.8) 

This form coincides with that obtained in [21-24]. Our derivation which 
has used (C.3) is more direct than prosented in [24]. Eq. (C.8) also coincides 
with the form originally suggested in [21] up to a change irrif — > {im — Aq — 
V)oo. This overall factor can be taken out to redefine a normalization of the 
partition function and one comes to (3.1). 

One can also show that nonzero modes can be neglected since the ratio 
{im — A„ — V)~^/{im — Aq — V)qq is of the order 0(mp). Namely for the 
two-channel situation, n = 0, 1 one has 

f \ m-i im -Xi-Vii - y -^ 2im-Xo 

(^--A-K)oo=^^^^^, (.m-A-V^), =^-^^ (C.9) 

Now take into account that Ai ~ Vu ~ p~^, while m and Aq are much smaller 
(see Appendix B). 
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Appendix D 

Study of the pion inverse propagator, N{k, k') in (5.6) 

We first prove that N[k = k' = 0) = 0. From (5.6) one has 

N(0,0) = -Tr\^^ iM + ^^ M J^ M] = (D.l) 

2 iD + iM iD + iM -iD + iM 

11 1 -^ 

= -Tr\—. M . D] = 

2 iD + iM -iD + iM 

1 1 -^ 1 -- 

= -Tr{ . D + -. D) 

4« ^ -iD + iM iD + iM 

Consider the tranformation of inversion P of all coordinates, Xa — * —x^ , B 



-B,. 



From (D.l) one can see that under P transformation A^(0, 0) changes sign. 
Therefore, when A^(0, 0) is integrated over all gluonic fields and is a number, 
it should vanish 

<Ar(0,0)>B=0 (D.2) 

From the invariance properties of the vacuum with respect to the shift of 
coordinates, one can deduce that < N{k, k') >b written as 

< N(k, k') >B= I du dv!e'^''+'^'''' < -Tr\-r^ imiu)6{u - u')+ (D.3) 

^ ' ^ i 2 \d + iM ^ ' ^ ' "- > 

iD + im -iD + iM 



where 



has the property 



eN 
<x\m{u)\y>= ^^^^K{x,y,u), (D.4) 



< N{k, k') >B= {27r)^5{k + k')N{k). (D.5) 

Expanding now A^(A;) and again using invariance of the vacuum we get 

A^(A;) = k'^Noik). (D.6) 
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Finnaly we report in this Appendix the calculation of < N{k, k') >b in the 
limit p — !► 0. We have for the first term in (5.6) 

N^^'^ =<-tr j dx dy du (D.7) 

fN 
■e'^'^'>S{x, y, B)^^^W^{y - u)^y, u, x)^+(x - u)tD >b 

When p ^ 0, both x and y tend to u and one can replace S{x, y; B) — > 
S{u, u; B). Integration over D{x — u) and {y — u) yields 

Ar(i) = {2nY5{k + k')const (D.8) 

Analysis of the second term on the r.h.s. of (5.6) yields in the limit p ^ 

iV(2) = fdudw 6*'="+^'^'"' < itr S{w, u; B)M{0)S{u, w; B)M{0) >b (D.9) 

Finally we can rewrite the sum A^*^^^ + A^(2) as 

<N{k,k')>B=N^^^ +N^^'> = (D.IO) 

= {27r)^6{k + k')const + M^{0)7r5{k)) 
This justifies eq. (6.10) used in the text. 
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